Many-Body Vortex Effect on the Transverse Force Acting on a Moving Vortex 



Nobuhiko Hayashi 
Computer Center, Okayama University, Okayama 700-8530, Japan 
(February 1, 2008) 



The nondissipative transverse force acting on one mov- 
ing vortex under the influence of another vortex is discussed 
in fermionic superfluid systems, where the relative velocity 
between the vortices is finite. On the basis of detailed nu- 
merical solutions of the Bogoliubov-de Gennes equation, the 
Berry phase for an adiabatic motion of the vortex line is ex- 
amined for a two- vortex system. It is found that the detailed 
electronic structure of a vortex core can affect the transverse 
force, without abandoning the previous discussions of the ro- 
bust Magnus force on a single vortex. 

PACS number(s): 47.32.Cc, 74.60.Ec, 74.60.Ge 

Much interest has been devoted to the dynamical mo- 
tion and pinning effect of vortices existing in varioyis 
situations in nature such as type-II superconductorsEl, 
supecEuidity in '^He, and superfluid neutron star mat- 
ters.Du It is essential also from the applied physics point 
of view that we understand transport properties due to 
vortices in type-II superconductors, especially for high- 
Tc cuprates, to better control superconducting charac- 
teristics under magnetic fields. In addition, the recent 
progress of experiments on '^He in a rotating cylinder 
shed new light on the fundamental physics field where 
exciting analogies between the. -early Universe and the 
superfluid ^He were examinedofl For these systems, the 
dynamics of moving vortices is the focal point. 

In this letter, we consider a many-body effect between 
vortices on the dynamics of a vortex. We will propose 
that a novel force should act on one moving vortex due 
to surrounding another vortex when the relative velocity 
between the vortices is finite. 

Currently, in the case of fermionic superfluids and su- 
perconductors, the fundamental equation of motion for a 
moving vortex is highly controversial. In particular, the 
term of the nondissipative transverse force in the equa- 
tion of motion, i.e., the vortex velocity dependent part 
of the Magnus force (per unit length_p£.the vortex axis 
z), az X Vi, is under intensive debateflia Here, is the 
velocity of a vortex. 

Broadly speaking, in the single vortex case, at present 
there are two conflicting points ofjideaj on the theory of 
the transverse force. One theoryE3~Hl considers a mo- 
mentum transfer from the (vortex) system under con- 
sideration to a heat bath system by a pilaxation of the 
quasiparticles of the vortex bound stateaij inside a vertex 
core. In this theory, called the spectral flow theory,Ej the 
coefficient of the transverse force, a, essentially depends 
on the vortex bound states (i.e., the electronic structure 



inside a vortex core) in combination with the relaxation 
time r of the quasiparticles. Experimental results on •^He. 
in a rotating cylinder support the spectral flow theory,E3 
which is the ccucial point of the ^He simulation of the 
early Universe.Q In contrast, the other theoryEjO claims 
that the transverse force on a moving vortex is a robust 
quantity which never depends on such quasiparticle re- 
laxation inside a vortex core. This "robust" Magnus force 
does not depend on the details of the vortex bound states 
inside a vortes_core, but only on the superfluid density far 
from the coralZI as long as a single vortex is considered. 
There also exist experimental results supporting the ro- 
bust Magnus force theoryjij and a negative experimen- 
tal result for the spectral flow theory is obtained there. 
Thus, the transverse force issue is controversial now, and 
careful interpretation of the experimental results by gen- 
eralizing the single-vortex theories mentioned above to 
more complicated situations will be required to resolve 
it, as noted in ref. ^. To advance our knowledge of vortex 
dynamics, it is important to consider many-body effects 
between vortices whose relative velocity is nonzero, due 
to, for example, vortex pinning. 

The transverse force is also a key point in the under- 
standing of the HalL effect in the mixed state in type- 
II superconductorsjij especially for the mysterious sign 
reversal of the Hall resistivity (conductivity) obsecyed 
experimentally in several type-II superconductors. cJ A 
typical objection to the robust Magnus force theory has 
arisen from small Hall angles observed in most transport 
experiments, because the robust (or full large) Magnus 
force leads to an extremely large Hall angle and to no 
sign reversal, as long as the many-body effect is not con- 
sidered. However, in real superconductors, the vortex 
pinning cannot be avoided, and then it is expected that 
the relative velocity of vortices is generally nonzetn, lead- 
ing to many-body effects between vortices. Aou previ- 
ously proposed that the movement of vortex vacancies 
in the pinned vortex lattice as a many-body vortex ef- 
fect should lead to a relevant result and the sign rever- 
sal, while he maintained the robust Magnus force acting 
on indis/idual moving vortices. Recent transport experi- 
mentsE3 certainly revealed that the vortex pinning played 
an important role in the Hall effect in some way. There- 
fore, it is required to further investigate what type of 
many-body effect can exist in various possible vortex sit- 
uations, (e.g., in the plastic flow state, in which portions 
of the jzprtex lattice move while other portions remain 
pinnedo), i.e., situations in which the relative velocity 
of vortices is generally nonzero. 

In addition, in high-Tc cuprate superconductors, it was 
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universally observed that the sign-reversal Hall effect ex- 
istedpOrily in underdoped samples, but not in overdoped 
ones.cil A hypothetical concept arises from this result; 
the change of the electronic state due to the doping could 
be related to the internal electronic structure inside vor- 
tex cores so that it affects the dynamic property of vor- 
tices. Such an idea opens a nea& possibility of relating the 
doping-dependent Hall effected to the electronic struc- 
ture of various vortex cores currently the focus of atten- 
tion, e.g., cores with the antiferroKiagnetic moment£j in 
dx2-y2 + idpy superconductivity,c3 and with doping de- 
pendence. EH 

Motivated by the above concepts, we investigate a con- 
tribution from the electronic structure around a vortex 
core to the transverse force acting on a moving vortex un- 
der the influence of another vortex (e.g., pinned one). We 
consider such a many-body effect between vortices whose 
relative velocity is nonzero, by discussing a two-vortex 
system. A Berry phase picked up by the system for an 
adiabatic vortex motion is examined to derive the coef- 
ficient of the transverse force az x v^, accorcliH|&4:a-.the 
Berry phase approach to the vortex dynamics ,ll3'c3'c3 on 
the basis of numerical wave functions around each vor- 
tex. It should be noted that the relaxation time r is 
not included in the present analysis, and therefore the 
contribution of the vortex core whicb-|We propose is in- 
dependent of the spectral flow theoryllj which essentially 
depends on t. We will consider a neutral system and rec- 
tilinear vortices (or vortices in a two-dimensional system) 
with vorticity antiparallel to z. 

We will basexLur analysis on the Bogoliubov-de Gennes 
(BdG) theory,E3 which is the spatially inhomogeneous 
version of the BCS theory. The system is described in 



terms of the Bogoliubov wave function (uj{r), 



r)). 



We start with the BdG equation given, in a dimension- 
less form, by 



K 
A*(r) 



Undir) exp[i{l - ^)9] and Vj{r) = Vn,i{r) exTp[i{l -\- ^)9] 
with A(r) = A(r)exp[ — 16*] in cylindrical coordinates, 
where n is the radial quantum number and the angu- 
lar momentum |^| = ■ ■. We expancLthe eigen- 
functions in terms of the Bessel functions O Jmi^) as 



Vn.l 



where 0i,„(r) = [V^/RJm+iia 
the z-th zero of Jm{r), 
R = 20-40^0- Equation 
trix eigenvalue problem, 
eq. (^, developed bj 
lized in many cases 
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and i = 1,2, • • • ,7V. We set 
|l]) is reduced to a 2A^ x 2A^ ma- 
This useful teclmique to solve 
l^nd Schliiter,E3 has been uti- 



Adiabatic vortex motion produces a Berry phase 0o in 



'j]Xj, and the total 
p system, |$) — > 



the solution of eq. (||), Xj ~^ exp[z 
Berry phase F picked up by the 
exp[zF]|$), is given as F = - J^j 

We consider the single vortex case first. When the 
center of the vortex ro{t) adiabatically moves, the mo- 
tion gives rise to the Berry phase (f>j picked up by the 
Bogoliubov wave function Xj ■ With the cylindrical sym- 
i n ^t|jr| of the vortex, the total Berry phase is obtained 

T = - J dt J d^ x{roVr J) S{r) 
= -7r[S'(oo) - 5'(0)] dt{xoyQ ~ yoio), (2) 



where S{r) is the canonical angular momentum density 
composed of Xji ^^id Tq = v/,. In the framework of 
the BdG theory, the physical quantity S{r) must be ex- 
pressed in the form with the Fermi distribution function 
f{E), as pointed out by GaitanJ^ namely 

S{r) = 2 [-f{E,u)\unAr)\'{-l + l 



E„j>0 



(1) 



+ {l-fiEn,l)}\VnAr)\'{-l 



(3) 



where K — —W^/2kp£^o — fi, /i is the chemical poten- 
tial, and ^o(=^^f/Ao) is the coherence length [Aq is the 
uniform gap at zero temperature T = 0, fcp (v-p) is 
the Fermi wave number (velocity), and h = I]. In eq. 
(|l|), the length (energy) is implicitly measured by £e| 
(Aq). The system is characterized by a parameter fcp^o-Eil 
The pair potential is self-consistently determined with 
A(r) = 5E|£,|<.o"jW";W{1 - Vim- Here, g is 
the coupling constant and wd the energy cutoff, which 
are related by the BCS relation via the transition tem- 
perature Tc and the gap Aq. We set wd — 20 Aq. 

To obtain exact solutions Xj(r) (i.e., electronic struc- 
ture) around a vortex, we numerically calculate eq. dJ) 
under the following conditions for clarity, as in ref. |3^. 
(a) The system is a cylinder with radius R. (b) The 
Fermi surface is cylindrical, (c) The pairing has isotropic 
s-wave symmetry. Thus, the system has cylindrical 
symmetry. We write the eigenf unctions as Uj{r) = 



At T = 0, eq. (^ which is calculated using our numer- 
ical solutions {uj,Vj) of eq. (|l|) gives S{0) = and the 
uniform value S — — /c/2 far from the vortex core (see 
Fig. |l|), where p/2 is half the total particle density and 
is, at T = 0, equal to the density of the Cooper pair, 
i.e., half the superfluid density Ps/2. Thus, the expres- 
sion for tiue Magnus force derived at T = by Ao and 
Thouless,E3 a = —hps/2, is obtained from eqs. (||) and 
(^ as in refs. ^ and |2^. In the present study, we extend 
the above discussion for T = to the finite temperature 
case within the mean-field BdG framework and with an 
ansatz that, at finite temperatures, the total Berry phase 
is the sum of the thermally weighted contribution from 
each Bogoliubov wave function, i.e., the total Berry phase 
should be obtained from eq. ^ by calculating the finite- 
temperature canonical angular momentum density S{r) 
composed of m^, vj, and f{Ej) (eq. (|3|)). It is consistent 
with an imaginary time path integral formulation.EJ We 
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should note here that the solutions of the BdG equation 
eq. (0), i.e., the Bogoliubov wave functions (uj,Vj), have 
all the information on the system on both the condensate 
and non-condensate in a two-fluid picture. 
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FIG. 1. The canonical angular momentum density S{r) 
and the total particle density p(r) around a vortex at sev- 
eral temperatures, in arbitrary units, but on a common scale, 
fcpfo = 4. Note the relation S = — p/2 at T = far from the 
vortex center r = 0. 



In Fig. |l|, we show S{r) at several temperatures. S(0) 
remains zero at finite temperatures in Fig. |l| due to the 
form of cq. (^) and the necessary structure of u„.;(r) and 
Vn,i{r) expanded by the Bessel functions. Because S{Q) — 
0, even at finite temperatures the transverse force does 
not depend on the vpttex core in the single vortex case. 
The particle densityP p{r) = '2Y.E,>o[\^ji'r)? fi.^]) + 
~ /(^i)}] ; is simultaneously presented in Fig. 
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rtj. Here, we adjust the chemical potential /i at each tem- 
perature so that the jiarticle density p far from the vor- 
tex core is invariableo With increasing temperature, \S\ 
decreases with respect to p/2. We can find from Fig. 
|l| that the decrease of \S\ far from the core obpvs the 
temperature dependence of the Yosida function.L3 This 
result gives numerical proof of the physical estimation in 
ref. 11^. Then, \S\ indicates the superfluid density Ps/2 
at finite temperatures. 

Let us consider the two-vortex problem, which is the 
main aim in the present study. We consider two vor- 
tices, labeled 1 and 2, moving with finite relative veloc- 
ity. Their relative dynamics is of interest here. Hence, 
without loss of generality, vortex 1 is at rest at position 
ri and vortex 2 at V2{t) moves around vortex 1. We as- 
sume the vortices are separated so that there is almost no 
overlap between the vortex bound states of each vortex. 
Then, wave functions bounded inside each vortex core 
are defined separately from each other. A two-vortex 
pair potential is expressed in the product form, A(r, t) — 
A„(r - ri)A„(r - r2{t)), where A„(r) = A(r)exp[-i0] 
denotes the pair potential of the single vortex. The wave 
functions near each vortex are expected to be close to 



those of the single vortex which have cylindrical sym- 
metry, but a phase factor is attached so that the above 
product form of the pair potential is constructed from 
those wave functions. 

Under the above conditions, w& abtain the bound-state 
contribution to the Berry phaseoO as F^ = rib+r2fc; the 
contribution from the bound-state wave functions of vor- 
tex 1 is Fih = ~2t:A J dt{R^Ry - RyR^)/\R\^, and that 
of vortex 2 is Fsb = -2t:A J dt[R^Ry - RyR^)/\R\^ - 
irSb f dt{x2y2 - y2i2), where R = r2{t) - ri, A = 
J rdrP{r), (the center of each vortex is r = 0), and 



P(r) = -2 J2 

0<E„j<A{T) 

+ \vnA{r)\'{l-f{E,,,i)} 



Unj{r)\^f{E„,i) 



(4) 



The second term in T2b is the bound-state contribution 
to the Berry phase which leads to the Magnus force in 
the absence of vortex 1, and is zero {Sb — 0) as seen 
in the single vortex case. What we aim to show in this 
letter is that the bound-state contribution P{r) is af^ 
tually nonzero. In the previous two-vortex discussion,c3 



this "lu," 



type" contribution of the bound states 



has been disregarded because of the equahty \uj\ = \vj\ 
for approximated bound-state solutions. However, ex- 
act solutions in general must have essential asymmetry 
between \uj\ and \vj\ due to the existence of the vor- 
tex itself, which leads to a local breaking of particle-hoia 
symmetry in the density of states inside the vortex coreEll 
and an electric charging of the core.E^ The Fermi func- 
tions should also be attached to the Bogoliubov wave 
functions on the analogy of S{r) in eq. In Fig. H, 
we show P{r) which is calculated with the exact solu- 
tions {uj,Vj) around the vortex, and find it is certainly 
nonzero. 
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FIG. 2. Plot of P(r) (in arbitrary units) obtained using 
exact wave functions of the vortex bound states around a 
vortex (see text), fcp^o = 4. 

On the other hand, the wave functions of the extended 
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states above the gap A(r) spread over vortices 1 and 
2. To consider thei^contribution, let us apply a discus- 
sion by DziarmagaE3 based on an analogy of the impu- 
rity .-effect on the Magnus force due to Ao and Thou- 
less.llj The Berry phase is proportional to the mean 
number of superfluid particles enclosed by a closed path 
around which a vortex is moved adiabaticallyEZI The sta- 
tionary vortex 1, as a kind of impurity, traps particles 
and excludes them from the extended states relevant 
to the Berry phases picked up by the extended wave 
functions. A quantity which becomes a central issue is 
6Se{r) — Se{r) + the deviation of Se from its 

asymptotic value —ps/2, where Se{r) is the contribution 
from the extended states to S{r) in eq. (||) around the 
stationary vortex 1. Then 2 J dr2TTr6Se{r) (= B) corre- 
sponds to the excluded particle number which works to 
decrease the magnitude of the|-S4ngle- vortex Magnus force 
I — hps/2\ in spatial average.E3 The contribution to the 
Berry phase is ST ^ -B J dt{R^Ry - i?2,i?^)/|R|2. 

Taking the variation of the Berry phase Ff, -f 5T about 
r2, we finally obtain the force F2 acting on vortex 2 due 
to the presence of vortex 1 , F2 = (/i(4yl+_B/7r)/|Rp) [zx 
R + z • (R X R)R/|Rp] , and we propose to add it to the 
general formEiJ of the usual vortex equation of motion. In 
addition to the nondissipative transverse force (the first 
term in large brackets), we obtain an interesting dissipa- 
tive term (the second one) which, to our knowledge, has 
not been reported so far. An analysis of this dissipative 
term is a future task. In what follows, we discuss the 
nondissipative transverse force. 
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FIG. 3. To the canonical angular momentum density S{r) 
around a vortex (in arbitrary units), the contribution from 
the vortex bound states Sb{r) and that from the extended 
states 5e(r), where S{r) — Sb{r) + Ss{r). The deviation 
SSeir) = 5e(r) + p,/2 (see text), fcp^o = 4. 

In Fig. 1^, we show our numerical results for Se{r) to- 
gether with the bound-state contribution Sb{r) and total 
S{r). The deviation dSe{r) is positive as noted in ref. |2|. 
It is of importance to note here that the deviation can 



become 6Se{r) > | — Ps/2|(= |S'(oo)|) locally around the 
vortex core, as seen in Fig.|3[ We are able to understand 
this by the factH that in general the extended states 
contribute in terms of negative rotation to the supercur- 
rent around a vortex, and the bound states contribute in 
terms of positive rotation to it. We note in Fig. ^ that 
at a higher temperature, the region of such large 6Se{r) 
expands following the divergence of the coherence length 
near T^. Accordingly, near Tc the region around such 
stationary vortices where dSe{r) > \ — Ps/2|, can pre- 
dominantly occupy the inside of, e.g., a superconductor, 
and thus the transverse force acting on the moving vor- 
tex can become not only a decreased Magnus force but 
also an opposite force. It may naturally explain the pSpn 
reversal Hall effect near Tc observed experimentally£3a 
The detailed structure of 6Se{r) can depend on the_de- 
tails of various vortex cores, and thus on the dopingpH in 
high-Tc cuprates. 

Effects of surrounding vortices on a vortex are usually 
taken into account through the superfluid velocity depen- 
dent part of the Magnus force, (3z x Vg, where is the 
superfluid velocity including the contribution of the sur- 
rounding vortices. In the present study, in contrast, we 
proposed quite a different force due to the surrounding 
vortices, az x , which acts on a moving vortex. Accord- 
ing to Newton's action-reaction principle, our transverse 
forces acting on the moving vortex simultaneously act on 
the surrounding vortices in each opposite direction. If the 
surrounding vortices are pinned ones, our force should be 
included in a depinning condition for those vortices. 

In conclusion, we considered how the details of the vor- 
tex core can affect the vortex dynamics by the many-body 
vortex effect other than the possibility of the spectral flow 
force which is due to r. The present study throws a light 
on the possibility for the transverse force which may de- 
pend on the details of vortex cores, without abandoning 
the robust Magnus force theory for the single vortex. 

A numerical simulation adopting the robust Magnus 
force and the force F2, and a direct two- vortex system 
analysis by numerically calculating the BdG for a system 
containing two vortices, would be interesting and are left 
for future work. 

I am grateful to Ping Ao for bringing my attention to 
his work and for helpful discussions. 
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